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Получены расчетные зависимости коэффициента теплоотдачи, скорости движения и трансцендентное уравнение для толщины 
ламинарной пленки конденсата, справедливые для регулярного теплового режима. 


В [1-3] рассмотрены задачи, связанные с расче- 
том стационарной теплоотдачи при конденсации 
пара с использования расчетных зависимостей ко- 
эффициента теплоотдачи, впервые полученные 
Нуссельтом. В условиях управления или регулиро- 
вания, например, отборов пара в отдельных ступе- 
нях турбины тепловой электрической станции про- 
цесс конденсации имеет нестационарный характер. 
Поэтому для такого процесса представляет практи- 
ческий интерес оценка нестационарной теплоотда- 
чи при пленочной конденсации пара на стенке. 


\Ѵ х (0,у) = 0, (4) 

д) граничные условия 

Т(т,у = о) = Т, (5) 

Т(т,8 х ) = Т, (6) 

ѴЛ(т,у = 0) = 0, (7) 

®ІА) =0 . (8) 
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Постановка задачи 


Пусть в процессе пленочной конденсации вся 
теплота, выделяющаяся на внешней границе плен- 
ки, отводится к поверхности охлаждения. В на- 
чальный момент времени движение пленки на 
стенке отсутствует, а вдали от стенки, т.е. на рассто- 
янии у = 8 Х (рис. [2], с. 354) изменение скорости не 
происходит. Перенос теплоты через пленку осуще- 
ствляется путем теплопроводности. 

Известна температура стенки, которая поддер- 
живается постоянной во времени, и она меньше по 
величине температуры насыщения - 7] при данном 
давлении. Принимается также известное допуще- 
ние [1] о том, что температура частиц на поверхнос- 
ти пленки конденсата равна температуре насыще- 
ния. Теплофизические свойства конденсата и пара 
считаются известными и постоянными величинами. 


Система уравнений, описывающая нестацио- 
нарный одномерный по координате процесс кон- 
денсации пара имеет вид: 

а) дифференциальное уравнение энергии 


Э Т _ Э 2 Т 
Эт “ а Э/ ’ 


( 1 ) 


т > 0; 0< у < <5 


б) уравнение движения несжимаемой вязкой жид- 
кости -уравнение Навье-Стокса 


ЭѴС 

Эт 


8 Х +Ѵ 


эх 


( 2 ) 


при следующих краевых условиях: 
в) начальные условия 

Т(т = 0, у) = Т 0 , (3) 


Здесь использованы известные обозначения [ 1 -3] . 

Решение задачи теплопроводности (1), (3), (5), 
(6), следуя [4-7], имеет вид 

Т(т,у) = Т +(Т - Т)-^-+— ^Г(«,т)8іпд у, (9) 
<5, тс „ = і 


где г (л, т) = - 

П 


і [г 0 [і-(-і)"] + 
[+[(—!)'" Г —Г] 


■ехр(-а/Тт), р п = пк / <5 . 


Аналогично находится решение системы урав- 
нений (2), (4), (7), (8), которое запишем в виде 


ѴѴ(т,>') = ^Х 

ѵ 


ло Х~ 

2^1 

2, — ехр (-ѵу„т ) 8іп у т у 

8 х „ = і у я 


, (Ю) 


здесь У,„ = (2 т - Г )к 1 2 <5 . 

Анализ решений (9), (10). Температура и ско- 
рость движения пленки конденсата являются 
функциями, зависящими от координат и времени. 
Стадия теплового регулярного режима наступает 
при числе Фурье Р 0 = ах/8\ > 0,25. Это означает, 
что в решении (9) можно пренебречь всеми члена- 
ми ряда за исключением первого. Пусть макси- 
мальная толщина ламинарной пленки конденсата 
8 Х = 1 ■ 1 0 -3 м. Температура насыщения Г, =127 °С и 
физические свойства воды [3]: Я = 0,686 Вт/(м-К), 
р ж = 939 кг/м 3 , а = 17, МО -8 м 2 /с. Тогда Р = 0,171т, 
т.е. процесс выравнивания температуры от 7] до Т с 
будет проходить в течение 1,5 с. На основе приме- 
ра [3] можно констатировать, что время наступле- 
ния регулярного теплового режима будет наблю- 
даться при т>5,2- 10^ с. Иными словами, нестацио- 
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нарный процесс конденсации пара является быст- 
ропротекающим процессом (табл. 1). 

Таблица 1. Изменение во времени т толщины пленки кон- 
денсата 5 Х и локального коэффициента теплоот- 
дачи сс х при х = 3 м 


т,с 

МО' 3 

5 -10' 3 

МО' 2 

5 -ІО' 2 

Стационарный 

режим 

5,, ІО" 5 М 

1,86954 

4,17691 

5,89155 

12,0618 

13,1307 

а _ Вт/(м'К) 

12,6 

138 

388 

3535 

5224 


Для стадии теплового режима изменение тем- 
пературы конденсата во времени будет подчинять- 
ся следующей зависимости: 

Т(т,у) = Т+ (Т-Т)—~ 

8 х 

2 ап 2 к у /,,4 

(Т + Т ) ехр(- — — г)8іп— . (11) 

п 8 8 х 

Тогда плотность теплового потока, согласно за- 
кону теплопроводности Фурье, будет равна 

Т - Т - 

ап 1 ■ (12) 
-2 (Т +Г)ехр(- — г) 

Если привести аналогию с явлением теплопро- 
водности [4] , то движение пленки конденсата, сог- 
ласно решению (10), происходит за счет увеличения 
во времени скорости движения пленки. При этом 
на твердой стенке ее скорость во времени равна ну- 
лю, а на поверхности, граничащей с паром, она не 
превосходит скорости невозмущенного потока. 

При значениях угл’ 2 /(4(5()>0,25 в решении мож- 
но пренебречь всеми членами ряда, кроме первого. 
При этом возникает погрешность, расчет которой 
оценивается по методу, изложенному в [4]. Для вы- 
ше рассмотренного примера с использованием 
данных [3]: ѵ = 0,24-1(Им 2 /с получим г = 0,45 с, т.е. 
наблюдается быстропротекающий процесс. Для 
регулярного режима средняя скорость движения 
пленки конденсата будет равна 

ІѴ( г) = — |У (г, ѵ)дѵ = ^-ЧѴ( г); 

* 8 2 х ѵ 

х О ѵ 


д = -А 


( дт \ 


Уду 


Я 

Л=о 8 , 


массовый расход конденсата через поперечное се- 
чение пленки шириной в 1 м запишется так: 

С(т) = р ж ІѴ х (т)8 х = 8 ' Р “ 8 ' 7У(т), 

V 

1 32 ( п ' 

здесь 2Ѵ(т ) = ехр ѵг . 

3 л \ 4<5; 

Если приравнять количество теплоты, выделяе- 
мое при конденсации пара, к теплоте, которая пе- 
реносится теплопроводностью к твердой поверх- 
ности стенки, можно получить уравнение для оп- 
ределения толщины пленки. Оно имеет вид 

Рж8 г8;Щт) = ХАТѵМ(т), (13) 


т +т 

где АТ = Т - Т , М(т) = 1-2 -*■ ^ 

5 с АТ 

Это уравнение решается методом последова- 
тельных приближений. 

Следуя [2], находим искомый коэффициент 
теплоотдачи: 

а = — М( г). (14) 

' 8 Х 

При стационарном тепловом режиме М(т)=1 и 
а х = Я / 8 х , 


\ 
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2 Л 
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здесь 


ЗЛА Тѵх 



Р ж ё/ 


Обсуждение результатов 

В качестве примера рассмотрим задачу 8-1 [3]. 

На поверхности вертикальной трубы высотой 
Н= 3 м происходит пленочная конденсация сухого 
насыщенного водяного пара. Давление пара 
Р = 2,5- ІО 5 Па. Температура поверхности трубы 
Т =123 °С. Необходимо определить толщину плен- 
ки конденсата 5 Х и значение местного коэффици- 
ента теплоотдачи а х в зависимости от расстояния х 
от верхнего конца трубы. При расчете следует счи- 
тать режим пленки конденсата ламинарным по 
всей высоте трубы. 


Таблица 2. Стационарная теплоотдача при конденсации пара на поверхности вертикальной трубы 


Координата х, м 

0,1 

0,2 

0,4 

0,6 

1,0 

1,5 

2,0 

3,0 

Расчет 
по [3] 


0,06 

0,0715 

0,0845 

0,094 

0,107 

0,118 

0,127 

0,140 

а х , Вт/(м 2 К) 

11430 

9620 

8150 

7320 

6530 

5880 

5410 

4900 

Расчет по формуле (13) 
м 

0,056105 

0,06672 

0,079346 

0,08781 

0,099772 

0,11042 

0,11865 

0,131307 

Расчет по формуле (14) 
а х , Вт/(м 2 К) 

12227 

10281 

8645 

7812 

6875 

6213 

5782 

5224 

Погрешность, % 

7,0 

7,0 

6,1 

6,7 

5,3 

5,7 

6,98 

6,6 
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В табл. 2 приведено сравнение 8 Х , а 1: . величины 
которых рассчитаны по приближенным формулам 
Нуссельта [1-3] и полученным в работе зависимос- 
тям (13) и (14). 

Из таблицы видно, что теория Нуссельта дает за- 
ниженные значения коэффициентов теплоотдачи по 
сравнению с соответствующими данными по пред- 
лагаемой зависимости (14). Отклонение не превы- 
шает 7 %, что можно объяснить заменой истинных 
значений на среднеинтегральные величины. Следует 
обратить внимание на "кажущуюся" высокую точ- 
ность проведенного расчета, особенно для нестацио- 
нарного процесса. Так, например, при х = 0,1 ми 
т = 1 • 1 0 3 с получены точные значения 

8 Х = 1,8677-10 5 м, а х = 360,4 Вт/(м 2 -К) при невязки 
между правой и левой частью уравнения (13) 
А = 1,84-Ю' 11 . 
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При 8 Х ~1,868. ІО- 5 м, а х ~ 304 Вт/(м 2 -К) 
(А = —2- 10 10 ), соответственно; если 8 Х = 1 ,87- 10 5 м, 
а х ~ -74,6 Вт/(м 2 К) (А = 8, 1- 10 10 ), что противоре- 
чит физическому смыслу. Поэтому все расчеты бы- 
ли приведены при невязке А < 2- 10 40 . Результаты 
расчетов показали, что характер изменения толщи- 
ны пленки и коэффициент теплоотдачи в стадии 
регулярного режима по высоте вертикальной трубы 
ничем не отличается от стационарного режима. 

Вывод 

Нестационарный процесс конденсации водяного 
пара при давлении менее 2,5 бар протекает в преде- 
лах от 1 до 5 мс. Показано, что в стадии регулярного 
теплового режима толщина пленки и коэффициент 
теплоотдачи возрастают в пределах 2-3 порядков. 
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